Boundary contributions to specific heat and susceptibility in the spin-1/2 XXZ chain 
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Exact low-temperature asymptotic behavior of boundary contribution to specific heat and sus- 
ceptibility in the one-dimensional spin-1/2 XXZ model with exchange anisotropy l/2<A<lis 
analytically obtained using the Abelian bosonization method. The boundary spin susceptibility is 
divergent in the low-temperature limit. This singular behavior is caused by the first-order contri- 
bution of a bulk leading irrelevant operator to boundary free energy. The result is confirmed by 
numerical simulations of finite-size systems. The anomalous boundary contributions in the spin 
isotropic case are universal. 
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Irrelevant operators in the renormalization-group (RG) 
theory are considered to give only subleading contribu- 
tions to any physical quantity by definition. A subtle ex- 
ception to this general rule can be found, however, when 
a bulk irrelevant operator has not so large boundary scal- 
ing dimension. Consider, for example, a Gaussian theory 
in (l+l)-dimcnsion with a bulk perturbing operator O 
having scaling dimension, x , larger than 2. First, let us 
impose periodic boundary condition to a one-dimensional 
(ID) system of length L. In second-order perturbation 
the irrelevant operator is found to give a contribution of 
order £jT 2Xo ~ 3 to the specific heat of the periodic system 
at temperature T . This is small compared with the lead- 
ing term of order LT. Next, let us suppose that the sys- 
tem has open boundaries. In general the specific heat has 
boundary contribution of order L . Indeed, first-order 
perturbation leads to a boundary contribution of order 
L°T X °~ 2 to the specific heat, where x' is boundary scal- 
ing dimension of the perturbing operator. If 2 < x' < 3, 
then the boundary specific heat gives a dominant con- 
tribution at low enough temperatures, even though the 
operator is irrelevant. 

This interesting possibility is in fact realized in the 
ID spin-1/2 XXZ model, as recently pointed out by 
Fujimotoi The model with open boundaries was stud- 
ied earlier by de Sa and Tsvelik using Bethe ansatz, 2 
who found anomalous boundary contributions to the spe- 
cific heat and spin susceptibility. Similar and more ex- 
tensive calculations 3 ' 4,5,6 for the spin isotropic case have 
shown that the boundary contribution to the spin sus- 
ceptibility at T = exhibits a Curie-like behavior with 
a logarithmic correction as a function of a weak exter- 
nal field. These results seemed somewhat mysterious, 
as irrelevant operators localized at a boundary can give 
at most a temperature-independent contribution to the 
susceptibility£2ai The origin of these anomalous bound- 
ary contributions was not understood, until Fujimoto and 
Eggert recently examined boundary contributions from a 
bulk leading irrelevant operator to the specific heat and 
the susceptibilityiifi Their theory seems rather compli- 
cated, however, as it uses properties of boundary states 
in the boundary conformal field theory. The aim of this 
paper is to present an alternative, simpler, and straight- 



forward derivation of the exact results using the Abelian 
bosonization method. We also verify our analytical re- 
sults with numerical calculations for the anisotropic XXZ 
model. 

Our results are relevant to many ID correlated elec- 
tron systems with impurities. An impurity potential in a 
ID system of repulsively interacting electrons is strongly 
renormalized and becomes perfectly reflecting in the low- 
energy limit^*iS*^ A long wire with a finite density of 
impurities is thus effectively cut into many pieces of fi- 
nite length. Measurements of the specific heat and the 
susceptibility of such a wire will exhibit the anomalously 
large boundary contributions. 

We analyze the spin-1/2 XXZ model in a weak mag- 
netic field, 
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where J > and L is the number of sites. We are inter- 
ested in the case where the exchange anisotropy A is in 
the parameter range 



1 



< A < 1, 
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for which we shall find an anomalous contribution to the 
boundary free energy. 

First we consider a semi-infinite spin chain, taking the 
limit L — > oo. Its free energy has, in addition to a bulk 
term, a boundary term of order 1 due to the presence of 
the boundary at j = 1. We calculate the boundary free 
energy at low temperatures using the low-energy effective 
theorj^ for the Hamiltonian QJ, which is the Gaussian 
model perturbed by a nonlinear operator. The Hamilto- 
nian density of the effective theory is 
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The lattice constant is unity. The bosonic field (f>(x) and 
its conjugate operator II(x), defined on the positive half 
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line, satisfy [<fi(x),H(y)] — iS(x — y). The parameter R 
and the velocity v are functions of A and h. At zero field 
h = 0, the parameter r\ defined by r\ — 2irR 2 is 
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with which the spin velocity v at h = is written, 
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We note that in the parameter range (J2J of our interest 
?y satisfies the inequality 



< V < 1 
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with the fields <fi and n. In particular, Sf is written as 
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where x — j and a is a constant calculated in Refs. fl5l 
andllfil For the semi- infinite wire the bosonic field obeys 
the boundary condition at x = 



<f>(x = 0) = 0. 



(8) 



Once Eqs. Q and J5J) are given, sign of the coupling con- 
stant A of the nonlinear operator is fixed to be positive. 
For the gapless regime (— 1 < A < 1) of the XXZ chain 
the exact value of A is known to b o 17 i 18 
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under the regularization condition on the zero- 
temperature correlator in the Gaussian model, 
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We are now ready to calculate the leading contribution 
to the boundary free energy. We begin with the path 
integral representation of the partition function Z, 
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where the Lagrangian density in the imaginary time is 



at h = 0. The original spin operator can be expressed £ 
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To eliminate the last term (oc h), we change the field as 
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' . 2nRv '. 
which yields a constant contribution —xoh /2 to the bulk 

free energy with the bulk susceptibility given by xo — 
L/(2irriv). This term is not of our main interest and will 
be discarded in the following analysis. 

Since the operator A cos(20/i?) is an irrelevant pertur- 
bation in the gapless phase, the low-temperature expan- 
sion of the free energy can be obtained by expanding the 
partition function in powers of A. The leading boundary 
contribution to the free energy then comes from the first- 
order term cx A in the expansion. After making the shift 
(jT3)l . the first-order term reads 
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Under the boundary condition JSJ the bosonic field can be expanded as 
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where uj n = 2irnT '. Inserting the mode expansion 1151) into Eq. (|14H and performing the Gaussian integral of <f>, we 
obtain 
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where Zo = Z\\ = q and we have introduced the high- (|10fl . The integral and summation in the exponent can be 
energy cutoff e~^ n ^ A with A = v to meet the condition easily performed, leading to the first-order contribution 
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to the free energy F = —TlnZ, which is nothing but the 
leading term in the boundary free energy, 
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where a is a short-distance cutoff of the order of half a lat- 
tice spacing. Note that F\ has a contribution only from 
the boundary region within the distance v/2itT from the 
edge of the spin chain. The exponent 2/rj is the boundary 
scaling dimension of the irrelevant operator cos(2</>/i?). 

Our next task is to evaluate the boundary free energy 
F\ at low temperatures. Let us first consider the simplest 
case h = 0. To obtain the low-temperature expansion of 
the integral in Eq. I|17[l . we change the variable as 



For a small but finite h we may expand cos(2hx /rjv) 
in Eq. (|T7|> as 
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where we have taken the limit a — > in the convergent 
integral. The integral in the right-hand side is elemen- 
tary. After some algebra we find 
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where ip'(x) = d 2 lnT(x)/dx 2 . The boundary contribu- 
tion to the susceptibility is then obtain as 
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where uq = tanh(27rTa/u) ss 2nTa/v. For the parame- 
ter range of our interest, the integral is divergent in 
the limit u Q — ► 0. We integrate I\ by parts to separate 
out the diverging contribution and obtain 
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where we have set uq = in the second term. The first 
term gives regular T 2n terms (n = 0, 1, 2, • • •) in the free 
energy, whereas the cutoff-independent second term leads 
to the anomalous contribution to the boundary free en- 
ergy. From Eqs. i|17[l and l|19|l we find the leading term 
in the boundary free energy 
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where E\ = Fi(T — 0) is the first-order boundary energy 
and subleading terms of order T 2 are ignored. Then the 
boundary contribution to the entropy is 
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and the boundary contribution to the specific heat is 
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For the XXZ chain with 1/2 < A < 1, Eqs. (g21> and 
@ give the exact leading boundary contribution to the 
specific heat. For A < 1/2 the leading term in the bound- 
ary specific heat is proportional to T. At A = 1/2 there 
appears a logarithmic correction to the T-linear specific 
heat. 



Xb-- 



d 2 F l 



dh 2 



h=0 



XT^na-l)^ 2 



2^'(i)] 



4^(2-^(2-1) 



— ) ■ ^ 



This is the exact leading contribution to the boundary 
susceptibility in the XXZ chain with 1/2 < A < 1 . Note 
that x is diverging as T — > 0. For A < 1/2 the leading 
term in the susceptibility is independent of temperature. 
A logarithmic correction shows up at A = 1/2. 
At T = the boundary energy E\ becomes 
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from which we obtain the boundary contribution to the 
total magnetization 
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where higher-order terms 0(h) are ignored. This is the 
exact form of the leading term in the XXZ chain for 1/2 < 
A < 1. The leading terms in C b (T), Xb(T), and M b (h) 
we calculated above are found in agreement with those 
obtained independently by Fujimoto and Eggert^ after 
some analytical transformations are made in their results. 

At the isotropic antiferromagnetic Heisenberg point 
(A = 1), we cannot substitute A of Eq. |J3J into Eqs. 
iPffil. (J2SJ), and lt2*7|) . as it vanishes at r\ = 1. To circum- 
vent this difficulty, we adapt a RG improved perturbation 
theory in which we replace a constant A with a running 
coupling constant X(E) at energy scale E. At h = 0, the 
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one-loop scaling equations for the sine-Gordon Hamilto- 
nian © are given by 



dl \ irR 2 
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where g = X/v and dl = —d\nE. The coupling g is 
marginally irrelevant at R = l/\/2~n (i.e., at A = 1) and 
decreases to zero logarithmically. To see this explicitly, 
we may write ^/2nR = 1 — e and expand the RG equa- 
tions H28f) in lowest order in e. We solve the simplified 
RG equations to find the solution corresponding to the 
antiferromagnetic Heisenberg chain, 
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for I ^S> 1. To obtain the leading term in C&(T), Xb(T), 
and Mf,(h), we may substitute 77 = 1 and 
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into the analytic formulas of the boundary free energy we 
obtained for 1/2 < A < 1. For example, the boundary 
entropy Sb is found from Eq. (|2~TJ> to be 
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from which the boundary specific heat becomes 
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in agreement with Ref. US Here T is a constant of order 
J. The boundary contribution to the susceptibility is 
obtained from Eq. I|25|l as 
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This result was confirmed later by the corrected calcula- 
tion in Ref. 0[ The boundary contribution to the mag- 
netization at zero temperature is found from Eq. (|27|l to 
be 
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in agreement with the exact result from the Bethe 
ansatz. 3-4 - 5 Here h$ is another constant of order J. We 
note that there exist subleading terms to these formu- 
las which are reduced only by a factor of l/ln(To/T) or 
l/\n{ho/h). Even though they are important for quanti- 
tative analysis, we do not try to evaluate it here, as it is 
beyond the scope of this paper. 

Finally, we present comparison between the ana- 
lytic result of Eq. (|26|l and numerical data of energy 




FIG. 1: E(L,m)/L as functions of 1/L for A = 0.8. The 
curves represent the fitting to Eq. 1351 . 
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FIG. 2: Numerical estimates of eo(m) (symbols) and their 
exact values (curves) as functions of m for several typical 
values of A. 



eigenvalues of the XXZ spin chains with open bound- 
aries. Using the density-matrix renormalization-group 
(DMRG) methodfi2*2S we have computed the lowest en- 
ergy E(L, m) of Ho in the subspace of the magnetization 
per site m = (l/L) 53j=i Sj. The maximum system size 
we calculated is L — 800. The number of block states 
kept in the DMRG calculation is up to 200. The numer- 
ical errors due to the DMRG truncation, estimated from 
the difference between the data with 150 and 200 states 
kept, are typically less than 10~ 5 . The numerical data 
are therefore accurate enough for the discussion below. 

The energy E(L, m) is known to depend on the system 
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FIG. 3: £i(m) as functions of (2-7rm) 2 for typical values of 
A. The curves are the fitting to Eq. I|360 . Inset: Numerical 
estimates of hi as a function of A (open circles) and the exact 
result, Eq. 13711 (solid curve). 

size L as 

E(L,m) = Le (m) + £ 1 (m) + Y^j. (35) 

The bulk energy £o(m) can be obtained exactly by solving 
Bethe ansatz integral equationSfSLSSiS while the bound- 
ary energy £i(m) is related to Eq. I|26l) . Here, we fit the 
numerical data to the formula (|35[) taking eq, E\, and £2 
as free parameters. The data of E(L,m) at A = 0.8 are 
shown in Fig. \I\foi several values of m. We have found 
that the fitting works pretty well for all As and ms calcu- 
lated. In fact, the estimated values of £0 obtained from 
the fitting are in excellent agreement with the exact val- 
ues calculated from the Bethe ansatz equations (see Fig. 

The bulk energy £o(m) has smooth quadratic depen- 
dence on 77i since the irrelevant operator cos(2-7r/i?) yields 
only subleading contributions to the bulk quantity. 

According to Eq. I|26|) and the relation between the 
bulk magnetization and the external field (m = xoh/L), 
the boundary energy £i(m) should depend on m as 

£l (m) = £i(0) - hi(2irm)$ -1 - (36) 

2X2 

for small m, where \i is a constant and the exact form 



of hi is 

AFQ-DsinKl-l)] 
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Note that the boundary energy of the open spin chains 
should be twice that of the semi-infinite wire that we 
have considered in the analytic calculation. To estimate 
the coefficient hi numerically, we fit the data of £i(m) to 
Eq. (|31)|) taking £i(0), hi, and % 2 as free parameters. We 
show the m dependence of £i(m) and the fitting curves 
in Fig. |3| Clearly, the plot of £i(m) versus (27rm) 2 shows 
large deviation from linear dependence on (27rm) 2 , con- 
firming the existence of the leading power-law term with 
the exponent 2/ry — 1, which emerges from the irrelevant 
operator Xcos(2<f>/R). The values of hi estimated from 
the fitting procedure are also shown in Fig. |3] We see 
that the numerical estimates agree extremely well with 
the analytic result (|37|l . This agreement can be regarded 
as a numerical proof of our analytic result, Eq. (|26|) . 

In summary, we have obtained the exact leading 
boundary contributions to the specific heat and the sus- 
ceptibility using the low-energy effective theory for the 
spin- 1/2 XXZ model. Our method can be easily extended 
to other ID systems (e.g., the extended Hubbard model 
and spin ladders in high magnetic fields) as long as their 
low-energy effective theory is given by Eq. J2J , which con- 
tains nonuniversal constants v and A to be determined 
with some other method, however. Our formulas 132L 
and f° r the spin isotropic case (A = 1) are 
universal, in the sense that if we make any moderate de- 
formation of the isotropic antiferromagnetic Heisenberg 
model that respects the SU(2) and translational symme- 
try and that docs not generate a spin gap, then these 
formulas apply to some values of To and ho, except for 
cases where the coupling A is fine tuned to vanish. This 
implies that, for example, the ID (extended) Hubbard 
model and other systems of interacting electrons with 
the spin SU(2) symmetry should universally exhibit the 
boundary contributions GUI, and 
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